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DETAILED ANALYSIS AND TEST CORRELATION OF A 
STIFFENED COMPOSITE WING PANEL 

D.D. Davis, Jr. 

Abstract 

State-of-the-art nonlinear finite element analysis techniques are evaluated by applying them to 
a realistic aircraft structural component. A wing panel from the V-22 tiltrotor aircraft (shown in 
Fig. 1) is chosen because it is a typical modem aircraft structural component for which there is 
experimental data for comparison of results. From blueprints and drawings supplied by Bell 
Helicopter Textron, Inc., a very detailed finite element model containing 2284 9-node Assumed 
Natural-Coordinate Strain (ANS) elements was generated. A novel solution strategy which 
accounts for geometric nonlinearity through the use of corotating element reference frames and 
nonlinear strain-displacement relations is used to analyze this detailed model. Results from linear 
analyses using the same finite element model are presented in order to illustrate the advantages and 
costs of the nonlinear analysis as compared with the more traditional linear analysis. Strain 
predictions from both the linear and nonlinear stress analyses are shown to compare well with 
experimental data up through the Design Ultimate Load (DUL) of the panel. However, due to the 
extreme nonlinear response of the panel, the linear analysis was not accurate at loads above the 
DUL. The nonlinear analysis more accurately predicted the strain at high values of applied load, 
and even predicted complicated nonlinear response characteristics, such as load reversals, near the 
observed failure load of the test panel. In order to understand the failure mechanism of the panel, 
buckling and first-ply failure analyses were performed. The predicted buckling load was 17% 
above the observed failure load while first-ply failure analyses indicated significant material 
damage at and below the observed failure load. 



Figure 1 . Wing Panel From V-22 Tiltrotor Aircraft 


ORIGINAL PAGE 

BLACK AND WHITE PHOTOGRAPH 


INTRODUCTION 


Complex structural systems which incorporate advanced materials and novel configurations 
require accurate and dependable structural analysis techniques. The finite element method has 
emerged over the years as the leading analytical tool of the engineer. However, the nonlinear 
response of very complex structural components has been difficult to predict. Recent advances 
have made such problems tractable. Advanced finite elements have been developed that are less 
sensitive to modeling distortion such as warp, taper, skew, and extreme values of aspect ratio than 
traditional element formulations (Refs. 1-3). Nonlinear solution algorithms have been developed 
that accurately predict response of highly nonlinear, post-buckled structures (Refs. 4-8). In 
addition, advancements in computer hardware and software make it possible to predict the 
response of large structural components whose finite element models require thousands of degrees- 
of-ffeedom (Ref. 9). 

A major activity is underway at NASA Langley Research Center in structural analysis methods 
research to provide better analysis tools and techniques for the aerospace industry. With the need 
for advanced structural analysis methods in mind, a computer code has been developed that serves 
as a framework for structural analysis methods research. This research code, referred to as the 
Computational MEchanics Testbed (COMET), allows new technology to be implemented and 
compared with other methods in a common software system (Refs. 10-12). COMET currently 
contains state-of-the-art technology such as improved finite elements, equation solvers, 
eigensolvers, and solution strategies. There is also ongoing research in constitutive modeling, 
error detection and control, failure methods, global/local analysis, adaptive mesh refinement, and 
multiple methods interfacing techniques. 

New analysis techniques are often verified using simple textbook problems that have well- 
known theoretical solutions (Ref. 13). However, even methods that perform well on simple test 
cases may not be feasible for analysis of realistic structures. The primary purpose of this paper is 
to examine the performance of recent developments and new techniques by applying them to a 
complex realistic structure. In particular, a novel finite element formulation and a nonlinear 
solution algorithm are assessed, and both are described in a later section. A more detailed version 
of this paper as well as results from mesh convergence and verification studies are presented in 
another publication '. 

APPROACH 


In order to demonstrate the analysis techniques in a manner that is meaningful to industry, a 
structure was chosen that is typical of a modem aircraft structural component. A panel from the 
lower wing skin of the V-22 tiltrotor aircraft (see Fig. 1) satisfies this requirement. The panel is 
made from graphite-epoxy composite material and contains design features such as ply drop-ofts, 
ply interleaves, axial stiffeners, transverse ribs, clips, brackets, and a large central elliptical access 

hole. 


Design blueprints and working drawings of the panel, as well as experimental data, were 
supplied by Bell Helicopter Textron, Inc. From the blueprints, a very detailed finite element model 
of the panel was generated. The entire panel, including the clips, brackets, stiffeners, nbs, and 
load-introduction fixtures, was modeled with shell finite elements. Offsets were used to model the 
eccentricities caused by the ply drop-offs and interleaves. Boundary conditions were used to 
simulate the test fixtures and the applied loading. The finite element model is discussed briefly m a 
later section and in detail in another publication >. 


tA NASA Technical Paper entitled "Detailed Analysis Of A Stiffened Com^site ^g Panel 
Finite Element Modeling, Analysis Techniques, and Test Correlation by W. J. Stroud, D. D. 
Davis, Jr.,T. Krishnamurthy, and S. L. McCleary which comprehensively describes the analysis 

project is currently under preparation. 
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Many simplified models of the panel, as well as different discretizations of the full model, 
were analyzed in this study. The results of these preliminary models, although not discussed in 
this paper! provided great confidence in the element formulation and the nonlinear solution 
algonthm. Based on results of the preliminary models, the benchmark finite element 
generated. Using this finite element model, linear and nonlinear stress analyses were performed. 
The analytical results are compared to the strain gage data. The linear and nonlinear results are also 
compared with each other to illustrate the costs and benefits of the nonlinear analysis. 

PANEL DESCRIPTION 


General Description 

The panel of interest is a large, 35 inches wide by 78 inches long, all-composite panel from 
the lower wine skin of the V-22 tiltrotor aircraft. Figure 2 is an illustration of the test panel 
showing the various components and regions of interest. There are five I-shaped stiffeners, t e 
center stiffener is discontinuous due to the presence of the elliptical access hole. The access 
hole is large enough (12 inches by 21 inches) to permit visual inspection the interior q{ the wmg. 
This panel is a major structural component having a design ultimate compression load of 334,000 
Ss and a design limit compression load of 260,000 pounds Because of the heavy engines 
mounted to the ends of the wing, the panel also has specified stiffness requirements. In order to 
meet the strength and stiffness requirements, several design features were incorporated that greatly 
complicate dialysis of the panel. The skin beneath the stiffeners is padded l up by m ^kaving 0^ 
plies into the basic skin. The region surrounding the access hole is also padded up, but with ±45 
olies instead of 0° plies. The ±45° plies help to transfer the load around the cutout to the adjacent 
stiffeners The relatively "soft" skin (90% ±45° plies) along the edge of the cutout minimizes the 
stress concentration at the free edge where delaminations of the composite material are likely to 
occur. A complex system of graphite-epoxy and metallic test fixtures are attached to thepanel at 
each end of the access hole to simulate the bulkhead -type transverse nbs of the actual wingbox. 



Figure 2. Sketch of Key Components of Wing Panel 
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The Skin 

The panel skin is very complex due to changes in the stacking sequences caused by numerous 
ply drop-offs. Basically, the skin can be described in terms of five distinct regions. Figure 3a is a 
schematic of the panel with the ribs and stiffeners removed and illustrates the five regions of the 
panel skin. Figures 3b, and 3c are schematics that illustrate the ply stacking sequences for cross- 
sections AC and AB, respectively, from Figure 3a. 

Region 1 is termed the basic skin and consists of 19 plies. Region 2 is the padded region 
beneath the stiffeners, and is merely the basic skin with 27 0° plies strategically interleaved into it. 
Figure 4 is a schematic showing a cross-section of the padded region of the basic skin. Region 3, 
which surrounds the access hole, contains 46 plies and is essentially the basic skin with 27 ±45° 
plies interleaved into it. The layups and thicknesses of Regions 1 , 2, and 3 are constant over the 
regions and are explained in detail in another publication^. However, Regions 4 and 5 are 
transition regions, i.e., regions of the skin where the layup changes point to point due to insertion 
or replacement of plies. Region 4 is the transition region that transforms Region 1 into Region 3 
by interleaving 27 ±45° plies. Thus, the thickness of Region 4 gradually increases from 19 plies 
(at point A) to 46 plies (at point B), as is seen in Fig. 3c. Region 5 is the transition region that 
transforms Region 2 into Region 3 by replacing 27 of the 0° plies of Region 2 with ±45° plies. 
Thus the thickness of Region 5 remains constant (46 plies) while the layup changes from 
predominantly 0° plies at point C to predominantly ±45° plies at point D. A more detailed 
description of the stacking sequences and layups for Regions 4 and 5 is found in another 
publicationt. 


POINT D 
POINT C 
POINT A 


POINT B 


REGION 1 BASIC SKIN 

g 1 I region 2 PADDED SKIN BENEATH STIFFENERS 
REGION 3 THICK SKIN AROUND HOLE 
■ lit rim REGION 4 TRANSITION region of basic skin 
[TJTJTJTJ REGION 5 TRANSITION REGION OF PADDED SKIN 

a) 5 Sections of Skin 




b) Padded Region Beneath Stiffeners c) Ramp Transition Region 

Figure 3. Description of Skin Components 


t Ibid 
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The Stiffeners and Ribs 

The panel contains five identical stiffeners, four of which (Stiffeners 1, 2, 4, and 5 of Fig. 2) 
extend the full length of the panel. The center stiffener (Stiffener 3 of Fig. 2), interrupted by the 
elliptical access hole in the center of the panel, is divided into two 20-inch sections. The construc- 
tion of the stiffeners consists of five graphite/epoxy components: two identical back-to-back C- 
channels (15 plies thick), a cap (14 plies thick), and two filler strips. The filler strips are merely 0° 
plies rolled up to form a plug for the void areas caused by the radius of the comers at the top and 
bottom of the junction of the back-to-back C-channels. 

The caps of the stiffeners are relatively stiff, with approximately 60% 0° plies. The 
dimensions, stacking sequences, and reference directions for all the components of the stiffeners 
are given in another publication^. The entire stiffener is cocured and then cobonded to the skin 
with a precured strip of cross-ply cloth between the skin and stiffener. Stiffeners 1, 2, 4, and 5 are 
cobonded to the top of the padded regions of the basic skin, while Stiffener 3 is cobonded to the 
top of the transition region (Region 5 of Fig. 3a) along the centerline of the panel. 

The bulkhead-type transverse ribs of the actual wing panel are simulated in the test panel by a 
combination of L-shaped graphite-epoxy channels and metallic test fixtures which are attached to 
the panel at each end of the access hole. This system is bolted to the caps of Stiffeners 1, 2, 4, and 
5, attached to the reinforced skin between Stiffeners 2 and 4 by a graphite-epoxy shear tie, and 
attached to the end of Stiffener 3 by two L-shaped graphite-epoxy clips. The purpose of the 
simulated ribs is to provide out-of-plane support for the panel at each end of the discontinuous 
center stiffener, as do the bulkhead-type ribs of the actual wingbox. Detailed diagrams with 
dimensions and stacking sequences for the rib system are presented in another publication '. 

TEST DESCRIPTION 

Testing of the panel was performed by the Lockheed Georgia Company. The test panel was 
equipped with 84 uni-axial strain gages and 8 rosettes, for a total requirement of 108 data 
acquisition channels. The gages were placed in strategic locations over the entire panel, including 
both sides of the skin, the caps of the stiffeners, and along the edge of the access hole. Figure 4 is 
a photograph of the panel in the testing machine. The strain gages and strain gage wires are seen 
attached to the exterior side (skin) of the panel. The stiffener (interior) side of the panel was 
equipped with a comparable number of strain gages. 

The instrumented test panel was aligned in the testing machine, and a compression load was 
applied to the test panel by lowering the upper crosshead. The horizontal struts shown in Fig. 4 
provided out-of-plane support at each end of the access hole. Pins and metallic angle brackets 
were used to attach the struts to the rib systems. The panel was loaded to failure (405 kips), with 
strain gage data recorded every 52 kips up to limit load (260 kips), and then every 26 kips until 
failure. Strain gage locations and the resulting test data were provided to the authors in the form of 
diagrams, tables, charts, and graphs. 

THE ANALYSIS TOOLS 


The Software Framework 

The analyses were performed using the Computational MEchanics Testbed (COMET). 
COMET is an advanced structural analysis code developed at NASA Langley Research Center to 
provide a framework for structural analysis methods research (Refs. 11-12). This research code 
contains modules called processors, which are independent FORTRAN executables that perform 
specific functions. Equation solvers, eigensolvers, mesh generators, element stiffness matrix 


t Ibid 
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processors, system matrix assemblers, and post-processing packages are all examples of 
processors. COMET also contains a high-level command language (CLAMP) that, along with the 
modularity of the processors, provides the user with complete control over the analysis, thus 
enabling analysis methods research. CLAMP contains many FORTRAN-like functions such as 



Figure 4. Panel in Testing Machine 


logicals, looping, variables, and math functions. By using CLAMP to control the input to and the 
execution of the processors, very complicated parametric studies can be performed. New 
techniques can be developed and assessed in the COMET environment without developing all of 
the supporting software. For example, new equation solvers or finite element formulations can be 
put into COMET as processors and assessed while using the remainder of the software system and 
utilities. This allows researchers to concentrate on their areas of expertise without spending large 
amounts of time generating redundant supporting software. 

The Finite Element Formulation 

Although there are 10 families of finite elements in COMET, the finite element formulation 
used in this study is a 9-node assumed natural-coordinate strain (ANS) shell element, denoted 
within COMET as EX97. The EX97 element has five degrees-of-freedom per node; there is no 
"drilling" degree-of-freedom. The element is "shear-flexible", i.e., the formulation includes 
transverse shear stiffness terms. The formulation is based on the assumed strain in the natural 
(element) coordinate system. A detailed description of the element formulation is given in Refs. 1- 
3. As compared with other, more traditional formulations, the EX97 proved to be superior, 
requiring fewer elements for convergence and giving more accurate results when the mesh was 
distorted 1\ The nonlinear strain-displacement relations of the EX97, when used in conjunction 
with the corotational capabilities (Ref. 10) of COMET, have given very good results, even in 
analyses that involved extremely nonlinear, large-rotation, structural response. 


t Ibid 
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Solution Strategies 

t inpar Stress Analysis The results for the linear stress analysis were generated by solving the 
traditional linear system of equations, F=KU, where F is the load vector of external applied 
forces K is the assembled linear stiffness matrix of the finite element system, and U is the vector 
of unknown displacements. Based on the solution for the unknown displacements, the stress 
resultants and strains are then calculated at the element Gauss integration pointS T The stress 
resultants and strains are extrapolated to the nodes using the element shape function. In the results 
section of this paper, the strain results are reported at the nodes, since the location of the certain 
nodes were made to coincide with the placement of strain gages on the test panel. 

ripnmpirirallv Nonlinear Stress Analysis Although COMET contains many variations of nonlinear 
solution techniques, the results from only one such technique are presented herein. TJe nonlinear 
solution strategy of interest is a corotational Newton-Raphson algorithm with linearized 
Crisfield/Riks arc-length control (Ref. 8). The linearized equations of motion are solved iteratively 
until based on the convergence criterion, the converged solution is found. Since the tangent 
stiffness matrix is updated only at the beginning of each new load step, the technique is referred o 
as a modified Newton-Raphson method. An approach associated with continuation methods for 
nonlinear problems is based on controlling an equilibnum-path-arc-length parameter. Arc-lengt 
control techniques have been developed primarily by Riks (Refs. 14 and 15), Wempner (Ref. 7), 
and Crisfield (Refs. 4 and 5). The nonlinear solution strategy, the convergence criterion, and the 
arc-length control strategy are described in detail in Ref. 8. 

Bifurcation Buckling The buckling results reported herein were generated by solving the 
traditional linear structural stability problem, K4>; + A.;K G CTj - 0. where K is the 

assembled linear elastic stiffness matrix of the system. The assembled geometric stiffness matrix, 
Kr dependent only on the state of stress, can be calculated based on the linear or the nonlinear 
deformed shape of the structure. The i-th eigenvector, <fc, of the equation/solution is the i-th mode 
shape of the buckled structure. The i-th eigenvalue, X/, is the i-th buckling load factor, t.e., t e 
multiple of the applied load that would cause bifurcation of the solution. 

First-Plv Failure. Analysis For the frst-ply failure analyses, a maximum strain failure criterion was 
used with the material allowables listed in Table 1. In a first-ply failure analysis, a linear or 
nonlinear stress analysis is first performed. The state of stress is then calculated for each layer o 
everv element. The layerwise strain components are then compared to the material allowables, it 
any of the material allowables, by definition the ply has failed. Thus, a first-ply failure analysis is 
a "snap-shot" analysis, and is based on the state of stress at a given applied load. When a ply 
"fails" the loads in the ply are not redistributed to adjacent plies and the stiffness ot the 
ply/laminate is not reduced as is done in progressive failure techniques. Because of this, in the 
literal sense, first-ply failure analysis is only valid up to the point when the first ply fails (hence the 
term "first-ply failure"). However, with an understanding of these limitations, a rough (though 
nonconservative) estimate of the damage to the composite material at a given value of structural 

load can be made. 


MATERIAL 

STIFFNESS (MSI) 

Vl2 

FAILURE STRENGTH 

S(n in/ln 

) 

Eli 

E22 

G12 

Xt 

Xc 

Yt 

Yc 

s 

T300/5208 

21.0 

1.35 

0.95 

.34 

8290 

8290 

3880 

23880 

9480 

CROSS-PLY 1 
PI OTH 

3.27 

3.27 

5.40 

.34 

3880 

23880 

3080 

23880 

9480 

v- Lv< 1 1 1 

ALUMINUM 

io.o 

10.0 

3.54 

.30 

X 

X 

X 

X 

X 


* Failure was not considered for the aluminum fixtures. 
Table 1. Material Properties Table 
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FINITE ELEMENT MODEL 


Although many simplified preliminary models and perturbations of the benchmark finite ele- 
ment model were studied, only the results for the benchmark model are presented herein. The 
entire finite element model, close-up views of the edge of the cutout and the edge of a transverse 
rib, and schematics of the models of the transition region and stiffeners are shown in Fig. 5. 
Several element types were also studied, but only results for the most robust formulation, the 9- 
node ANS element, are presented herein. Based on the results of many preliminary models and 
mesh convergence studies, the relatively coarse mesh around the cutout was determined to be 
adequate. Another publication t contains more detailed descriptions of the preliminary models and 
results of the convergence studies, the modeling sensitivity studies, and the element performance 
studies. 

The skin, longitudinal stiffeners, and transverse ribs are modeled with the 9-node ANS 
elements discussed previously. The padded regions and the regions where the stiffener flanges 
were bonded to the skin were modeled as one contiguous composite finite element through the 
thickness. Therefore, the appropriate offsets were specified for the elements in order to maintain 
the eccentricity of the structure. Figure 5 shows a schematic of the stiffener and the corresponding 
finite element model. The bold lines represent edges of shell elements; the black dots represent the 
nodes. The eccentricity of the structure is modeled by specifying offsets for the appropriate 
elements. 

The benchmark finite element model contains 2284 elements (all 9-noded ANS elements), 
9486 nodes, and 4928 degrees-of-freedom. There are a total of 46 different material section 
properties. Since the objective was to predict the experimental results, the analysis was performed 
for only one load case and one constraint case, as is described in the next section. 

Composite Material Properties 

Due to the complexity of the composite material layups that comprise the panel, only a brief 
description is included in this section. More details about the stacking sequences for the individual 
components are given in another publicationt , where the exact material layups of the test panel are 
illustrated and the corresponding section properties of the finite element model are described. 

In COMET, each section property is defined by specifying the number of layers in the section 
and the material, the angle of orientation, and the thickness for each layer. A laminate analysis 
utility within COMET (processor LAU) then performs the through-the-thickness integration to 
calculate the smeared orthotropic laminate properties (the A, B, and D-matrices of Ref. 16). The 
laminate matrices include the effects of offsets in the shell walls arising from the built-up structure. 

The majority of the finite element model of the panel is comprised of multi-layered T300/5208 
composite material; the thickness of each layer is .005 inches. A single-layer strip of precured 
cross-ply cloth was cobonded between the stiffener and the padded skin. It was modeled as one 
layer in the contiguous layup containing the padded skin, the strip, and the flange of the stiffener. 
Some parts of the transverse ribs, though actually metallic test fixtures, were modeled with finite 
elements having the properties of steel. The layerwise material properties are given in Table 1. 

The ply drops of the transition region are modeled with four section properties as opposed to 
modeling all 27 longitudinal ply drops with an individual section property. The layups used for 
these four sections are the layups from the actual panel at the midpoint of the individual sections. 
A schematic of the ramp region (from point A to point B of Fig. 3), and the corresponding finite 
element model are shown in Fig. 5. The appropriate section thicknesses and offsets are also 
indicated. 


t Ibid 
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CROSS SECTION Of STIFFENER 



Figure 5. The Benchmark Finite Element Model 
Loads and Boundary Conditions 

To simulate the support struts of the test fixtures, the out-of-plane displacements are 
constrained to be zero at the midpoint of the top of the fixture brackets (see Fig. 5). At these 
points, rotations about the global X-axis are constrained and rotations about the global Y-axis are 
free, to simulate the pins that attach the struts to the rib fixtures. Since the test article is potted at 
both ends, providing a virtual clamp, all degrees-of-freedom are constrained to be zero at the end 
(X ^0) of the panel. To simulate the compressive motion of the testing machine crosshead, 
uniform end- shortening displacement is specified at X = L, while all other displacements at X = L 
are constrained to be zero, as is shown in Fig. 5. The "drilling" degree-of- freedom (the rotational 
degree-of-freedom whose vector is normal to the plane of the element) is constrained for every 
appropriate node, since the finite element formulation does not have these "drilling" degrees of 
freedom. 


ANALYSIS RESULTS & TEST CORRELATION 
Global Response 

Displacements The load-shortening curves for the linear and nonlinear analyses, as well as the 
failure load of the test panel, are shown in Figure 6. The global end-shortening response, even in 
the nonlinear analysis, is very linear through failure. Hence, the panel maintains its stiffness even 
when heavily-loaded. The nonlinear effects caused by the eccentricity of the discontinuous center 
stiffener have negligible effect on the global stiffness of the panel, at least up to the failure load. 
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The deformed geometry plots from the linear and nonlinear analyses are shown in Figs. 7a 
and 7b, respectively. The deformations are scaled by the same relative amount to allow qualitative 
comparisons to be made. Although the linear and nonlinear analyses indicate approximately the 
same deformations in the ribs, the out-of-plane deformation at the edge of the cutout is significantly 
greater in the nonlinear analysis. 

As the compressive forces develop in the panel, the skin attempts to expand in the transverse 
direction due to the Poisson effect. However, the caps of the stiffeners areattached to the relatively 
stiff transverse ribs, which prevent this Poisson expansion. The result is a complex state of stress 
in the region between the two transverse ribs. Due to the eccentricity caused by the discontinuous 
center stiffener, the region around the cutout between the two transverse ribs bends upward, 
causing an out-of-plane deflection in the positive Z-direction. Because of the complex construction 
of the test panel, the deformation pattern between ribs is m=l (one half-wave in the X-direction) 
and n=2 (two half-waves in the Y-direction) across the panel. Thus, the two stiffeners along the 
edge of the panel (Stiffeners 1 and 5) are bending in the opposite direction from the two stiffeners 
adjacent the access hole (Stiffeners 2 and 4). 



Figure 6. Load Versus End- Shortening 



a) Linear Analysis b) Nonlinear Analysis 

Figure 7. Deformed Geometry Plots 
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Local Response 

Strains At Panel Midlength Figure 8 is a plot of the axial strain (£\) across the panel midlength 
for an applied load of 156 kips. The strain gage data are shown as well as strains on the bottom 
surface of the skin and the top surface of the caps of the stiffeners from the linear and nonlinear 
analyses. There were 6 strain gages placed along the midlength of the bottom (skin side) of the 
panel, one at the centerline below each of the Stiffeners 1, 2, 4, and 5, and one at each edge of the 
cutout, as is indicated in Fig. 8. Four strain gages were placed along the midlength on the stiffener 
side of the panel, one each at the midlength centerline of the top surface of the caps of Stiffeners 1, 
2 4, and 5. At this value of applied load, both the linear and nonlinear analyses correlated well 
with test data. Notice the high gradients of strain near the edge of the cutout and across the width 
of the stiffeners. The high gradients in the stiffeners is an indication of bending in the plane of the 
panel, as can be seen pictorally in Fig. 7. 
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Figure 8. Axial Strain Across Panel Midlength, P=156 kips 

Figures 9a and 9b are plots of the axial strain (£ x ) across the panel midlength for an a PPjied 
load of 396 kips, which is very near the observed failure load of the test panel. Figures 9a and 9b 
show values for the bottom surface of the skin and the top surface of the caps of the stiffeners from 
the linear and nonlinear analyses, respectively. At this value of applied load, the linear analysis 
(Fig 9a) did not show good correlation with the strain gage data. The linear analysis predicted that 
the strain on the cap of Stiffeners 2 and 4 would be slightly higher than the strain in Sttffem 1 
and 5. However, the strain gage data indicate that the strain in Stiffeners 2 and 4 is about half of 
the strain in Stiffeners 1 and 5 at an applied load of 396 kips. The nonlinear results of Fig. 9b 
much more accurately predict the strain at this load. 
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Figure 9. Axial Strain Across Panel Midlength For Top and Bottom Surfaces of Skin 
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Point Strains as F„ no tion of Applied Load The axial strain (Ex) is plotted as a function of applied 
loadinFig 10. Figure 10a shows results at the midlength of Stiffener 1 for Jeeap and the 
bottom of the skin beneath the stiffener. Similarly, Fig. 10b is for Stiffener 2, and Fig. 10c is for 
*e top and bottom surfaces of the skin at the edge of the hole. The open circles represent discrete 




c) Top and Bottom of Skin at Edge of Cutout 

Figure 10) Axial Strain as Function of Applied Load 

lnad cteos from the nonlinear analysis, the dashed lines represent the extension of the linear path, 
antf i^fin^^mbu'irrepresent strain gage data. The vertical l* n f at 4^5 kipsindi^t^^i^to^d^at 
whirh the test oanel failed The divergence in the strain between the top and bottom surfaces ot an 
etmento smcurewUh increasing applied load is an indication. of bending. Stiffener 1 
conhn^stobendas is indicated by the initial divergence of the strains in Ftg. 10a, up toaloadof 
4(Xnips At 400 kips, the stiffener actually begins to unbend, or straighten, as is indicated by the 
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converging strains of Fig. 10a. If the panel had not failed at 405 kips, the nonlinear analysis 
predicts that, at a load of 453 kips, the stiffener actually would have began bending in the opposite 
direction. At the same load where Stiffener 1 begins to straighten (400 kips), the bending of 
Stiffener 2 becomes severe, as can be seen pictorally in Fig. 7b and numerically in Fig. 10b. 
Stiffener 2 is adjacent to the cutout and is the primary mechanism by which load is redistributed 
around the cutout. Initially Stiffeners 2 and 4 are bending in opposite directions than Stiffeners 1 
and 5, however when the bending of Stiffeners 2 and 4 becomes severe, it actually begins to 
straighten Stiffeners 1 and 5. The nonlinear analysis predicts that, had the panel not failed at 405 
kips, by the time the applied load is 453 kips, all four load-canying stiffeners would be bending in 
the same direction. Although the thick skin beneath the stiffeners carry most of the load, the axial 
strains at the edge of the cutout are the highest in the panel, as is seen in Fig. 10c. These high 
strains at the edge of the cutout play an important role in the ultimate failure of the panel. 

In general, there is excellent agreement between the strain gage data and the stress analysis 
results. The linear analysis accurately predicted the strains through the DUL of the panel. The 
nonlinear analysis results slightly under-predict the amount of bending in Stiffeners 1 and 2 at high 
loads, but the important response characteristics discussed above are correctly identified. The 
strain at the edge of the cutout, the most critical region in the panel, was accurately predicted all the 
way through failure of the panel. Although the linear analysis correlated well at the lower loads, 
the nonlinear analysis was necessary to predict the complicated response of the heavily-loaded 
panel as it approached the failure load. 

EXPLANATION OF FAILURE 


Buckling Analysis 

As was described previously, linear and nonlinear buckling analyses were performed. The 
predicted linear buckling load is 495 kips, approximately 17% higher than the observed failure load 
of the test panel. Contrary to what was expected, the nonlinear buckling load (the buckling load 
calculated by considering the nonlinear deformed geometry in the geometric stiffness matrix K c ) 
was not significantly less than the linear buckling load. The critical buckling mode is shown in 
Fig. 11. Notice that the buckling mode is a simple mode, with all stiffeners bending in the positive 
Z-direction. 



Figure 1 1. Critical Buckling Mode 
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First-Ply Failure Analysis 

Linear and nonlinear first-ply failure analyses were performed as was described previously 
In Fig. 12a and 12b, respectively, the elements are shaded according to the percentage of failed 
plies for the linear and nonlinear analyses for an applied load of 415 kips (the test panel failed at 
405 kips). Notice that at this load, both the linear and nonlinear results indicate a significant 
number of elements with plies that have strain levels exceeding their allowables. The nonlinear 
analysis predicts more elements with damage, but the linear analysis predicts a much higher 
percentage of failed plies per element. These differences are attributable to the fact that the 
nonlinear analysis predicts more bending in the region around the cutout. As the material is 
compressed and bends, the axial strain on one surface is intensified by the bending and the strain 
on the other surface is relieved. Therefore, the nonlinear analysis predicts damage at a lower load, 
or in this case more widespread damage for a given load, and the linear analysis over-predicts the 
damage as the entire thickness of the element is compressed beyond the compression allowables. 



a) Linear Analysis b) Nonlinear Analysis 

Figure 12. Percentage of Failed Plies at 415 kips Applied Load 

Discussion of Failure Analysis Results 

The first-ply failure results indicate that failure strain allowables were exceeded in a number of 
plies in the region around the cutout at the observed failure load of the test panel. Because of this 
and since the linear and nonlinear buckling loads were about 17% above the failure load, it is felt 
that the panel failed due to excessive strains in the region around the ctuout. The nonlinear analysis 
predicted the most significant damage to be at the edge of the cutout, a slight angle off of the 
vertical. The photograph of Fig. 1 3 shows that the failure zone passed through the edge of the 
cutout a slight angle off of the vertical. Due to the catastrophic nature of the failure, it is not known 
that the damage initiated at the edge of the cutout. These conclusions are preliminary at this point 
because sensitivity studies have not yet been performed to determine the sensitivity of the results to 
the material failure strain allowables that were used. It is also unknown at this time how sensitive 
the buckling and nonlinear results are to the imperfections of the stiffeners, because the model 
assumes that the stiffeners are perfectly straight along the length. 

SUMMARY AND CONCLUSIONS 

New finite element analysis techniques were evaluated by applying them to a complicated 
composite wing panel from the V-22 tiltrotor aircraft. A detailed finite element model with a 
relatively coarse mesh of 9-node elements was generated, and linear and nonlinear stress analyses 
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were performed. At low values of applied load (up to the design ultimate load of the panel), the 
linear stress analysis accurately predicted the strains and structural response characteristics of the 
panel. However, a nonlinear analysis was required to accurately predict the very complicated 
nonlinear structural response as the applied load approached the failure load of the panel. 
Calculated surface strains were very accurate as compared with strain gage data, even in critical 
regions of the panel such as the edge of the cutout and the caps of the stiffeners. Since the first-ply 
failure analysis indicated damage at the edge of the cutout at the observed failure load of the panel 
while the calculated buckling load was 17% higher than the failure load, it appears that material 
damage caused by excessive strain initiated the failure of the panel. 



Figure 13. Photograph of Failed Test Panel in Region Around Cutout 
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